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© Raman spectroscopy of respiratory gases. 

© Method and apparatus for determining the com- 
position and concentration of gases (14) present in a 
patient's airway by measurement of the spectrum of 
Raman scattered light (16) from these gases. The 
gases (14) present are assumed to be drawn from a 
predetermined set of gases with known Raman scat- 
tering spectra, and the concentrations are deter- 
mined by solution of a matrix equation Ac = b, 
where the c vector components are the unknown 
concentrations and the b vector components are 
determined from measurements of the Raman scat- 
tering intensities in a plurality of wavelength or 
wavenumber intervals. The linear system of equa- 
tions represented by the matrix equation may be 
overdetermined and may require solution by non- 
conventional methods, such as singular value de- 
composition. 




BNSDOCID: <EP 0557658A1_I_> 



Rank Xerox (UK) Business Services 

G. 10/3.5/3.3. i» 



r 

1 EP 0 557 658 A1 2 



Field of the Invention 

This invention relates to determining the com- 
position and concentration of an arbitrary mixture of 
gases in a patient's airway, and more particularly to 
such determination by means of Raman spectro- 
scopy. 

Background of the Invention 
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Measurement of the composition of a gas mix- 
ture is especially important in the airway of an 
individual who is connected to an auxilliary breath- 
ing apparatus. Examples include ventilators which 
assist respiration in the intensive care unit and 75 
coronary care unit in hospitals, and anesthesia ma- 
chines in hospital operating rooms ("ORs"). These 
patients are often monitored closely for their vital 
signs, including respiratory gas exchange. 

In the operating room, the anesthesiologist 20 
supports the patient's respiration and also controls 
the patient's depth of anesthesia using special gas- 
eous agents. A patient airway will contain the nor- 
mal respiratory gases, such as oxygen (O2), carbon 
dioxide (C0 2 ), and water vapor (H 2 0), or a mixture 25 
of nitrous oxide (N 2 0) and one or more haloge- 
nated anesthetic agents ("HA", usually halothane, 



separate measurement, using a paramagnetic or 
polarographic sensor, is needed to measure O2, 
which has no infrared ("IR") spectrum. Also, 
separate cells are needed to measure the HAs 
and the CO2 + N 2 0 concentrations. In addition, 
these devices have no means to detect other 
gases which may be present. Such devices may 
introduce error in the measurements, and the 
devices cannot be adapted easily to measure 
new agents. 

(2) Mass Spectrometry. Mass spectrometers 
("MS") can provide extremely accurate mea- 
surements of gas concentrations. Historically, 
MS devices are expensive and complicated in- 
struments that require frequent calibration and 
maintenance. These devices require use of a 
delicate vacuum system and ion source. Typi- 
cally, many operating rooms share a single MS. 
This reduces the response time between mea- 
surements, and requires relatively long sampling 
lines, which can distort the gas samples. In 
addition, the mass spectra of the various OR 
gases are not unique. Nitrous oxide and carbon 
dioxide have the same mass number, and 
isoflurane and enflurane are isomers. Therefore, 
one must observe fragmentation of these mol- 
ecules in the system and employ special al- 
jrithms to distinguish the molecules. The sys- 



the system. Various metabolic products and special 30 
gases, for example, to measure pulmonary func- 
tion, may also be present. Also, interferences from 
extraneously introduced gases such as ethanol and 
isopropanol may be present. Because severe injury 
can result from use of an improper gas mixture, 35 
anesthesiologists prefer to measure the composi- 
tion of the patient breathing mixture. Every compo- 
nent of the mixture is important. Measurement of 
oxygen concentration helps prevent hypoxia. Pres- 
ence of the CO2 waveform indicates healthy gas 40 
exchange. Measurement of the type and concentra- 
tion of the various anesthetic agents helps control 
and adjust the proper depth of anesthesia. The 
presence of other gases can indicate leaks and 
possible system malfunctions. 45 

Several Anesthesia Monitors ("AMs") exist to- 
day that perform this function. Four different tech- 
nologies compete for most of the market today. 
(1) Infrared Absorption plus Oxygen. One class 
of AMs uses infrared absorption to measure the 50 
halogenated anesthetics, C0 2 , and N 2 0. While 
this technique is widely used, it has the dis- 
advantages that it is difficult to distinguish the 
HAs from each other, because they have similar 
IR spectra that must be measured in the region 55 
X = 3-4 urn. Identification of the individual 
agents requires measurement in the far infrared 
region. X > 10 urn, which is more difficult. A 



gases, including helium and water vapor. Sharec 
systems cannot provide continuous, breath-to- 
breath analysis of airway concentrations, which 
is desirable. At least one manufacturer has intro- 
duced a small, stand-alone MS, offering each 
OR suite a dedicated measuring device. Al- 
though this device eliminates the problem of 
multiplexing measurements for different users, 
and has a more convenient design, the device 
still requires protection from contaminants, and 
can only be configured to measure a few spe- 
cific gases at a time. The problems of high cost 
and distinguishing between isomers and be- 
tween identical mass spectra remain. 

(3) Photo-acoustic Spectroscopy. These devices 
also utilize the properties of infrared absorption 
to characterize gas mixtures. A precision micro- 
phone detects pressure waves which are pro- 
duced when the gas sample absorbs IR energy. 
The sound level indicates the concentration. 
Like the IR systems described above, these 
devices need a separate system to measure 
oxygen concentration, and the technique cannot 
easily distinguish different HAs from each other. 

(4) Raman Spectroscopy. Scattering of light by 
the Raman effect has received much attention 
from scientists since its original exposition by 
C.V. Raman in 1928. Simply stated, when mon- 
ochromatic light illuminates a vibrating molecule, 
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light scatters in a process which decreases or 
increases the frequency of the scattered light by 
exactly the vibrational frequency of the mol- 
ecule. The shift in frequency of radiation is 
characteristic of the scattering medium, and is 
independent of the frequency of the illuminating 
radiation. Thus, measurement of the Raman- 
scattered light can be used to infer the prop- 
erties of the medium, such as the chemical 
composition and concentration. For measure- 
ments of OR and airway gases, this technique 
has the advantage that each OR gas, including 
oxygen and any poly-atomic molecule, has a 
unique Raman spectrum. 
Additionally, the Raman spectrum for a gas is 
usually contained in a relatively narrow wavelength 
band, which simplifies detection. Thus, Raman 
spectroscopy offers the promise of simultaneous 
measurement of all airway gases with a single 
measurement and less complex technology. 

In Raman scattering, a small fraction of colli- 
sions of photons with an atom or molecule are 
inelastic, with a photon either giving up a small 
portion of its initial energy E 0 to the collision part- 
ner and scattering as a photon of reduced energy 
E < E 0 (Stokes waves) or the collision partner 
giving up a small portion of its initial energy so that 
the photon scatters with increased energy E > E 0 - 
(anti-Stokes waves). In Rayleigh scattering of a 
photon with an atom or molecule, by contrast, the 
energy of the scattered photon is equal to the initial 
energy of the photon. This does not include light 
that is absorbed and re-emitted by processes such 
as phosphorescence or fluorescence. In a typical 
scattering situation, the ratio of intensity of 
Rayleigh to initial light intensity for gases might be 
about 1CT 9 and the ratio of intensity of Raman 
scattered light to initial light intensity might be 
about 1CT 12 . The change in wavelength for a Ra- 
man scattered photon of initial wavelength X 0 = c/ 
fo and scattered light wavelength X R = c/ f R is 
given by 

AX = x R -Xo = (c/f R )-(c/fo), (1) 

where AX > 0 (AX < 0) corresponds to appearance 
of Stokes waves (anti-Stokes waves). Substantially 
all scattered light at moderate initial energies arises 
from Rayleigh scattering or Raman scattering. For 
molecules of moderate or higher symmetry, not all 
modes of molecular vibration result in Raman lines. 
Some of the modes of molecular vibration produce 
infrared absorption lines but not Raman scattering 
lines, some modes produce Raman lines but not 
infrared lines, some modes produce both Raman 
and infrared lines, and some vibration modes pro- 
duce neither Raman nor infrared lines. 
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The Raman scattering cross-sections are ex- 
tremely small, and the intensity of the scattered 
light is very weak, as noted above, especially in 
gases where the molecular number density is also 

5 relatively small (compared to liquids and solids), 
and are therefore difficult to measure with accu- 
racy. The differential intensity of the Raman-scat- 
tered light scattered into a differential solid angle 
dn along a differential path length dz in a single 

w component gas is given by the formula 

dPfiaman = P 0 n 0 (60/6Q) dO dz. (2) 

where P 0 is the intensity of the incident light, n G is 
75 the number density of the scattering molecules, 
and (da/dQ) is the differential scattering cross-sec- 
tion in a given direction. The direction of the scat- 
tering is also dependent on the polarization of the 
incident light. If the gas contains more than one 
20 component and the components do not interact 
appreciably with one another, the intensity of each 
Raman line of a gas component is proportional to 
the concentration of that component so that Equa- 
tion (2) above can be used with a small modifica- 
25 tion to take account of the presence of the other 
components. 

For a given intensity of the incident radiation 
and sample concentration, one can maximize the 
magnitude of the measured Raman signal only by 
30 increasing the solid angle of the light-collecting 
optics, or by increasing the observation path length 
(i.e., using a larger scattering volume), because the 
molecular properties of the sample are not variable. 
In a conventional system for observing Raman- 
35 scattered light, a laser beam is brought to a focus 
in the medium of interest, creating a minute region 
of relatively intense electrical field, which excites 
the Raman effect. The light scattered from this 
region is collected by an optical system, typically a 
40 simple lens which images the scattering region 
onto a suitable optical filter and detector. The dif- 
ficulty of improving this simple design is evident 
from the observation that the etendue (defined as 
the product of the collection area and collection 
45 solid angle) is conserved in any ideal optical sys- 
tem. Thus, increasing the solid angle of light col- 
lection (lower f number optics) decreases the ob- 
servable area (and thus the path length). Each 
experimental system must optimize these param- 
50 eters against its own constraints. Some Raman 
systems further improve signal intensity by provid- 
ing multiple passes of the incident light beam 
through the observation volume, effectively increas- 
ing the path length, or by placing the observation 
55 volume inside an optically resonant cavity, effec- 
tively increasing the incident intensity P Q . 

Analysis of Raman scattered light is especially 
useful where the Raman spectrum of each of the 
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components present in a sample is a relatively 
simple line pattern and the Raman lines of the 
different components do not coincide with or lie 
close to one another. A complete range of vibra- 
tional frequencies can be covered with one moni- 
toring instrument, and the sample container can be 
glass or many other relatively transparent materi- 
als. Water may be present: the Raman spectrum of 
water is weak and diffuse in the band 200-3300 
cnT\ but the spectrum has a strong, broad peak 
centered at 3652 cm" 1 . The approximately linear 
relationship between component concentration and 
Raman-scattered light intensity makes the calcula- 
tion of concentration straightforward. Integration 
across a portion of the Raman spectrum to deter- 
mine the intensity of specified lines is also straight- 
forward. However, the sample should be non-flu- 
orescent, and the sample to be analyzed should be 
relatively transparent, with little or no absorption at 
the wavelengths of interest, and should be free of 
particulates. It is often difficult to apply Raman 
scattering to very low concentrations of the sample, 
because of the weak intensity of the individual 
Raman lines. 

Raman systems that measure the composition 
in patient airways must measure multiple spectral 
lines in order to distinguish ail the component 

ribe techniques 



disclosed by Tans et al in U.S. Patent No. 
4,630,923. Raman scattered light from the known 
and unknown gas samples is alternatingly received 
by a detector to determine the concentration ratios 
5 of two gases present in the unknown gas sample. 

A gas monitoring system, disclosed by Benner 
et al in U.S. Patent No. 4,648,714, collects light 
from a single focal region, illuminated by a properly 
oriented, polarized laser beam. A single set of 
w collection optics images the scattered light onto a 
single detector. An additional mirror, placed op- 
posite the collection optics, reflects Raman-scat- 
tered light back into the focal region and into the 
collection optics, further enhancing the signal. The 
75 sample cell can be inside the resonant cavity of the 
exciting laser to increase the signal still more. The 
system employs a rotating filter wheel that passes 
different Raman line filters in front of the detector 
so that different Raman lines are measured se- 
20 quentially. While this system can measure a larger 
number of Raman spectral lines (by adding more 
filters), the samples are not acquired simultaneous- 
ly from the same gas sample. This reduces the 
ability of the system to respond rapidly, as is 
25 desired for breath-by-breath analysis of the airway 
composition. 

In U.S. Patent No. 4,676,639, Van Wagenen 
discloses use of a gas ceil for Raman scattering 



Albrecht, in German Pat. DE 2723939C2, de- 
scribes a system of six detectors mounted in an 
equatorial plane around a region of focus. The focal 
region has a confocal cavity where the exciting 
laser beam makes multiple passes through the 
sample. This system requires six different detec- 
tion channels, each with associated collection op- 
tics, laser line rejection and Raman line filters, and 
separate detector. The configuration is limited to 
six channels and uses an unpolarized laser beam. 
More channels would be needed to measure all the 
OR gases of interest. In addition, the confocal 
cavity and each of the collection optics requires 
delicate adjustment to assure proper imaging of the 
focal region onto the detectors. 

In U.S. Patent No. 3,704,951, Chupp discloses 
use of a multi-pass gas cell for increasing the 
intensity level of light that is Raman scattered from 
a gas sample contained in the cell. Raman scat- 
tered light exits from the cell through a large side 
window in the cell. 

Leonard, in U.S. Patent No. 3,723,007, dis- 
closes a gas cell with a transparent side window for 
Raman scattering analysis and notes that two or 
more simple molecular gases may have distinct 
Raman shift spectra. 

A spectrophotometer that compares light Ram- 
nan scattered from a known gas sample with light 
Raman scattered from an unknown gas sample is 



30 The side windows may be coated wifhs^iarrow 
band anti-reflection coating, for passage of the Ra- 
man scattered light for detection outside the gas 
cell. 

Van Wagenen et al describe another system in 

35 U.S. Pat. No. 4,784,486, which uses multiple detec- 
tion units, each including collection optics (lens and 
back reflector), laser line rejection and Raman line 
filters, and a detector. Each detection unit collects 
Raman-scattered light simultaneously from sepa- 

40 rate focal regions. The system is similar to the 
system of Benner, but employs a completely sepa- 
rate detection unit for each Raman line, rather than 
inserting different filters serially. This design has 
the advantage that all the channels can acquire 

45 measurements simultaneously. However, the mea- 
surements are generated from different sample vol- 
umes and, thus, from different gas molecules. The 
sample gas must flow from location to location, and 
the flow rate must be such that all units see a 

so mixture of substantially similar composition. Thus, 
for a given flow rate and size of sample cell, the 
practical number of detector units is limited by the 
desired time response of each system. This cre- 
ates limitations because it is desirable to sample as 

55 small a volume as possible from the patient airway. 
An increase in the number of detection units in- 
creases the expense of the unit proportionally. 
Each unit also requires separate alignment and 
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adjustment, which increases the complexity and 
cost of the system. 

Many of the systems discussed above employ 
a combination of individual narrow-band filters. 
Each individual filter examines only a single spec- 
tral line or a small spectral region. This requires 
that a system use at least one filter for every 
molecule spectral peak of interest. Allowing for the 
expense of multiple channels and filters, these sys- 
tems are suited to the measurement of gas mix- 
tures, where the Raman spectrum of each gas 
consists of a single peak or a few well-separated 
peaks, and where these spectral peaks are unique 
for each gas in the mixture, with no overlapping 
spectral lines. In this case, each filter can measure 
each gas separately and completely. Overlapping 
spectral lines can be ignored if the gas of interest 
has another unique peak. This is true for the gases 
N2, O2, CO2, and N2O. The principal spectral lines 
of these gases are listed below as frequency shifts 
from the frequency of the illuminating source, such 
as a laser beam. Note that C0 2 and N 2 0 have an 
overlapping line at 1285 cm" 1 . 
N 2 , nitrogen : 2331 cm" 1 
O2, oxygen : 1555 cm" 1 
CO2, carbon dioxide : 1285 and 1388 cm" 1 
N 2 0, nitrous oxide : 1285 and 2224 cm" 1 
The limitations of these types of systems, 
which employ separate filters, arise when the Ra- 
man spectra of one molecule of interest becomes 
complicated, when it covers a larger spectral band 
or has many, scattered peaks, such as the haloge- 
nated anesthetics. These spectra are best char- 
acterized by their entire spectrum, not by single 
lines. The Raman spectra from the HAs have con- 
siderable spectral overlap with each other, and the 
broadband emission can contribute error to the 
signals of the respiratory gases and N 2 0. A repre- 
sentative mixture of the respiratory and anesthetic 
gases produces a complicated signal, with emis- 
sions and spectral lines across the entire spectrum, 
a mixture of all the individual spectra. Measure- 
ment of single spectral lines therefore gathers less 
signal and less information than a technique that 
gathers the entire spectral signature of the mixture 
or of a single component. Single filter systems use 
only small amounts of the total signal available, 
because such systems measure only part of the 
Raman emissions. For the HAs, measurement of a 
single line neglects most of the total signal. Single 
filter systems also neglect significant information 
about the Raman spectra. It is much easier to 
distinguish each HA from examination of its entire 
spectrum than from a single, possibly weak, line. 

A further limitation of the single filter system is 
that new or additional gases cannot be measured 
without the addition of new detection units with new 
filters. New gases which have substantial spectral 



overlap with existing gases may be very difficult to 
measure, even with additional filters. 

Also, if a Raman spectrum has closely spaced 
lines that must be distinguished from each other, 
5 the spectral bandwidth of the individual filters must 
be sufficiently narrow to differentiate nearby lines, 
or obscure undesired lines. This places a constraint 
on the spectral bandwidth of the laser source, 
which must be at least as narrow as the Raman 
10 line filters, and the laser must not drift in 
wavelength. This constraint also reduces the de- 
tected power. Further, filters can drift as they age. 
If a filter has a bandwidth of 1 nm, and the laser 
drifts in wavelength from 800 nm to 803 nm, the 
75 Raman line of interest moves out of the band of the 
filter and is not measurable. For these reasons, 
existing filter-based Raman systems often use gas 
lasers, such as argon-ion, which has narrow, con- 
stant lines and requires high input power. Solid- 
20 state lasers have wider line widths (> 1 nm), and 
often require special wavelength stabilization. 
Semiconductor lasers presently require low input 
power and produce wider linewidths. 

What is needed is a gas monitoring system 
25 that (1) allows simultaneous determination of the 
concentration of components in a multiple compo- 
nent gas; (2) is very efficient and uses as much of 
the total signal available as possible for examina- 
tion of the spectrum; (3) allows interrogation of a 
30 broad spectrum in a single measurement of a small 
single sample; (4) allows distinctions to be made 
relatively easily between gas components that may 
be present; (5) allows considerable freedom in the 
choice of sample container and sample concentra- 
35 tion; (6) allows the system to monitor and identify 
different compounds with minimal system modifica- 
tions; and (7) allows use of low power light sources 
and of relatively low cost detectors. 

40 Summary of the Invention 

These needs are met by the invention, which 
provides in one embodiment, a light beam of a 
selected wavelength \o, chosen to promote Raman 

45 scattering in at least one of a predetermined group 
of gases that may be present in the patient's 
airway; and a sample cell to receive and temporar- 
ily hold gas from a patient's airway, the cell being 
relatively transparent along the beam axis to light 

50 of wavelength \o and being positioned to receive 
the light beam so that gas in the sample cell is 
illuminated by the light beam. The invention also 
provides spectrograph means, having a beam en- 
trance and wavelength dispersion means, for re- 

55 ceiving a scattered light beam at the beam en- 
trance, for directing the scattered-and-received 
light beam to the wavelength dispersion means for 
dispersing this light beam into a plurality of light 
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rays, each ray having a different interval of light 
wavelengths. Light collection means is also pro- 
vided, for collecting portions of the light beam that 
undergoes Raman scattering by one or more of the 
gases in the sample cell and for directing these 
portions of the light beam to the spectrograph 
means beam entrance. The invention further pro- 
vides an array of photodetectors, each of which 
receives one of the plurality of wavelength-dis- 
persed light rays from the spectrograph means and 
issues a photodetector output signal representing 
the amount of Raman scattered light received by 
that photodetector. Finally, the invention provides 
computation means for receiving the photodetector 
output signal from each of M photodetectors (M > 
number of gases to be tested for) and for comput- 
ing the relative or absolute amount of Raman scat- 
tered light received by each such photodetector. 

This invention measures the complete Raman 
spectrum of a gas mixture to determine its com- 
position. The invention gathers information about all 
the gases present simultaneously, and from a sin- 
gle, small sample volume. Further, the invention 
will combine techniques of light collection, spectral 
separation, and signal processing to maximize the 
total amount of light collected by the spectrometer, 
and to utilize the information contained in the entire 



eters are often designed to observe a wide range 
of compounds, whose structure or spectra are often 
unknown. 

In contrast, this invention employs a grating 

5 spectrometer to measure the concentration, not 
molecular structure, of a limited number of known 
compounds, such as the airway gases. This inven- 
tion takes advantage of these differences to op- 
timize the combination of spectrometer optics, un- 

w conventional entrance slits, and diode array ele- 
ments for the determination of airway gas con- 
centrations. The optimized optical system, in fur- 
ther combination with an intensified source of 
Raman-scattered light, and an optimized data re- 

15 duction technique, are part of the present invention. 
Key elements of the new invention are use of 1) a 
means to maximize the intensity of Raman-scat- 
tered light from a small aperture, 2) a large eten- 
due for the spectrometer, and 3) an optimized data 

20 reduction means to determine concentrations. 

In order that the invention may be understood 
more easily exemplary embodiments will be de- 
scribed with reference to the following drawings in 
which: 

25 Figure 1 is a schematic view of one embodiment 
of the invention. 

Figures 2 and 3 are graphical views of the 
Raman scattering spectra for two gases, N2O 




of this invention to provide a means whereby new 
or additional gases can be detected in the airway 
gas mixture, without altering the basic spectrom- 
eter, sample cell or structure of the algorithm used. 

The invention is an efficient means for collect- 
ing the optical signals generated by Raman scatter- 
ing of an optical beam in gases. 
The basic invention consists of the following essen- 
tial components: 

1. A monochromatic source of illumination, e.g., 
a laser; 

2. A means of sampling and containing the 
gaseous media; 

3. A system of collection optics; 

4. A wavelength-selective filter to distinguish the 
Raman signals; 

5. A photodetector; and 

6. Computation means to determine gas com- 
position from photodetector signals. 

Earlier inventors have avoided using a grating 
spectrometer as a wavelength-selective filter be- 
cause of high cost and low optical efficiency. The 
invention described here combines several new 
ideas and technologies that reduce these limita- 
tions. Traditionally, grating spectrometers are used 
to observe the molecular structure of compounds. 
These systems are often optimized for high resolu- 
tion, and have low optical throughput - a narrow 
entrance slit and "slow" optics. These spectrom- 



patient's airway. 
Figures 4 and 5 are schematic views illustrating 
details of light illumination apparatus for the gas 
samples held in the sample cells. 
Figure 6 is a schematic view illustrating use of a 
35 variable wavelength filter to implement the in- 
vention. 

Figures 7, 8, 9 and 10 are schematic views 
illustrating use of light-emitting apertures with 
tubular sample cells according to other embodi- 
40 ments of the invention. 

A preferred embodiment 1 1 , shown in Figure 1 , 
includes an optical waveguide 13, through which a 
gas sample 14 flows. A laser beam 15 passes 
along the longitudinal axis of the waveguide 13 and 
45 produces Raman scattering along the entire 
waveguide length. The laser beam 15 is preferably 
directed through the waveguide 13 to minimize 
diffuse scattering of the laser beam from the inte- 
rior surfaces of the waveguide, but this is not 
50 essential. The laser that produces the laser beam 
15 may be a gas laser, such as He-Ne or argon- 
ion, but is preferably a solid-state laser diode emit- 
ting visible or near-infrared light in the wavelength 
range X = 500-850 nm. A preferred embodiment of 
55 the waveguide 13 is a cylindrical tube whose inter- 
nal surfaces are substantially fully reflective at the 
Raman shifted wavelengths of interest, for example 
a smooth silver-plated surface approximately 100 
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mm long with inside diameter between 1 and 2 
mm. In this manner, the effective length of inter- 
action (dz in Equation (2)) is increased substan- 
tially, thereby increasing the total Raman signal. In 
another preferred embodiment, the laser beam 5 
propagates in an optically resonant cavity that is 
surrounded by the waveguide 13. 

The scattered light is trapped in the waveguide 
or tube 13 and emits from each end of the 
waveguide in the same manner as that light emits 10 
from a solid optical fiber. The small diameter of the 
tube 13 provides a small aperture for imaging the 
Raman light 16 by the collection optics. The collec- 
tion optics preferably includes a lens 17 whose 
focal point is positioned near the aperture of the 75 
waveguide and that has an aperture with f number 
approximately equal to 1. The lens 17 substantially 
collimates the beam of Raman light and directs the 
beam through a laser line rejection filter 19 that 
blocks any scattered radiation at the wavelength of 20 
the laser beam 15. A second lens 21, also part of 
the collection optics, images the Raman light 16 
through the entrance slit 23 of a spectrometer 25. 
The spectrometer optics may also contain other 
elements that produce sharp slit patterns at the 25 
location of the spectrometer slit 23. The preferred 
embodiment has only one slit 23 at the entrance to 
the spectrometer 25. However, multiple slits, with 
different widths and spacings, may be used to 
increase the illumination inside the spectrometer 30 
25. 

The spectrometer 25 preferably includes a dif- 
fraction grating spectrometer 27. Raman-scattered 
light from the slit(s) 23 is directed at a concave 
diffraction grating 27, which is preferably a holo- 35 
graphic element or replica that is produced for high 
efficiency diffraction in the wave number range of 
interest (approximately 200 cm" 1 to 3300 cm" 1 
shift from the laser beam wave number). Use of 
replicas or masters of a holographic grating to 40 
disperse and image the Raman light reduces the 
system cost and system parts count, while provid- 
ing very good performance. The grating 27 dis- 
perses the Raman-shifted light 28, which contains 
all the Raman shift frequencies from all the spectral 45 
lines from the gas mixture, across the elements of 
a photosensor array 29. 

The photosensor array 29 is preferably a sili- 
con photodiode array of approximately 256 ele- 
ments, but can be a larger or smaller array of pin 50 
photodiodes, or a CCD array. Each element of the 
array 29 corresponds to a single, narrow 
wavelength band of light. Thus, the photosensor 
array 29 measures the entire spectrum of Raman 
scattering from the mixture simultaneously. The 55 
amount of light incident on each element of the 
array 29 is proportional to the amount of Raman 
scattering in that wavelength band and is received 



by computation means 30 that determines relative 
or absolute concentrations of the gases 1 4. 

The composition of the gas mixture can be 
determined by analyzing the measured spectrum, 
which represents the sum of the spectra of the 
individual gases, weighted by the concentration of 
each gaseous component. For illustrative purposes, 
the Raman scattering spectra for the gases N 2 0 
and C0 2 are shown in Figures 2 and 3, respec- 
tively, as a function of wavenumber shift A*' = VQ % - 
f R ' over the range 0 cm"" 1 S AS s 3,000 cm" 1 for 
Stokes waves. Each gas of interest may have one 
or several prominent spectral lines, and a promi- 
nent line of one gas may be adjacent to a promi- 
nent line of another gas. Other airway gases of 
interest include N 2 , 0 2 , SF 8 , halothane (CF 3 - 
CHCIBr), enflurane (CHCIF-CF 2 -0-CHF 2 ), isoflurane 
(CF 3 -CHCI-0-CHF 2 ), desflurane (CF3-CHF-O- 
CHF 2 ), sevoflurane (CH(CF 3 ) 2 0-CH 2 F), methoxyl- 
flurane (CH 3 -0-CF 2 -CHCI 2 ), ethanol(C 2 H 5 OH), 
isopropanol (C 3 H 7 OH), ether, cyclopropane, ac- 
etone (0 = C(CH 3 ) 2 ) and water vapor. 

Because the photosensor array 29 has sub- 
stantially more elements than the number of com- 
ponent gases, there exists a surplus of information 
from which to determine the concentrations: The 
solution is over-determined. Mathematically, one 
seeks the optimal solution vector c to the system 
of equations Ac = b, where A is a design matrix of 
the basis functions representing spectra of each 
individual gas, c is a vector representing the un- 
known concentrations of each gas in the mixture, 
and b is a vector representing the measured spec- 
trum from the diode array. The vector components 
b in the linear system of equations Ac = b serve 
as a forcing function whose values are determined 
by experimental measurement of the airway gas 
Raman scattering spectra. A preferred way to solve 
this problem and find the concentration compo- 
nents is to construct a Singular Value Decomposi- 
tion ("SVD") of the matrix A, and to determine a 
pseudo-inverse matrix A + , where c = A + b. If one 
uses this technique to determine the pseudo-in- 
verse matrix, then one solution is A + = (A ,r A) -1 A, 
where A ,r is the transpose of the matrix A, which is 
not necessarily square. Alternatively, one can con- 
struct a pseudo-inverse matrix A + using a least 
squares singular value decomposition ("SVD") ap- 
proach. 

Construction of a least squares SVD pseudo- 
inverse of a matrix A is discussed by G. Strang, 
Linear Algebra And Its Applications , Academic 
Press. 1980, pp. 137-145, which is incorporated by 
reference herein. For an m x n matrix A, the 
equation Ac = b may be an overdetermined sys- 
tem of equations because the number of spectra 
measurements (m) often exceeds the assumed 
number of distinct molecules (n), each of which is 
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assumed to correspond to a linearly independent 
spectral measurement in the vector space spanned 
by all such measurement values. The optimal 
choice of solution c = c minimizes the length ||Ac - 
b|| of the error vector, and one form of the solution 
is c = (A tr Ap 1 A tr b. This projects the forcing func- 
tion vector b into the sub-vector space spanned by 
the column vectors of the matrix A. The term 
Singular Value Decomposition refers to the de- 
composition of the matrix A into the form 



A = Qi EQ2, 



(3) 



where Qi and Q2 are m x m and n x n orthogonal 
matrices, respectively, and I is an m x n diagonal 
matrix with a sequence of positive diagonal entries, 
followed by a sequence of zero-valued diagonal 
entries. The diagonal entries of the matrix E are 
determined by solving an eigenvalue problem, ap- 
plied to A tr A, and the matrices Qi and Q2 are then 
determined by a Gram-Schmidt orthonormalization 
process. The SVD process is known, but its ap- 
plication to the spectral measurement versus con- 
centration problem appears new. 

Advantages of this embodiment include the 
following: 

1. All concentration of the components of the 

ale time and 



the center and focal point of the parabolic mirror. 
Two opposed ends of the sample cell 33 can be 
mirrors 35 and 37 that form part of an optically 
resonant cavity for the laser 39, as illustrated in 

5 Figure 4. 

The two opposed ends of the sample cell can 
also be mirrors that form part of a non-resonant 
cavity 41, where the laser beam makes multiple 
reflections through a small focal region, as illus- 

10 trated in Figure 5. These embodiments have all the 
advantages of the embodiment shown in Figure 4, 
except for the increased interaction length, and 
have a potentially narrow, linear scattering region 
that images exactly into the slit of the spectrom- 

75 eter. Also, this configuration may decrease the 
amount of unwanted scattered light into the spec- 
trometer. 

Any of the embodiments above may use mul- 
tiple slits at the entrance to the spectrometer. Mul- 
20 tiple slits cause multiple spectral images to appear 
on the diode array. The images are substantially 
identical but are translated on the array by a dis- 
tance determined by the spacing of the slits. The 
images sum together on the array to produce a 
25 new, mixed signal. For example, if a spectrum 
contains a single Raman line, two slits will create 
two line images, spaced apart, on the detector 
array. Because these "multi-slit spectra" are still 



2. The sample volume can be potentially quite 
small and is determined by the volumes of the 
waveguide, resonant cavity, or confocal cavity; 

3. The system has a small number of parts and 
no moving parts - one set of imaging optics, one 
diffraction grating, and one photodiode array; 

4. The aperture of the waveguide or cavity is 
small, and use of multiple slits increases the 
illumination in the spectrometer; 

5. Measurement of the entire spectrum of the 
mixture increases the total optical signal, pro- 
duces a better signal-to-noise ratio, and pro- 
vides extra information for more accurate cal- 
culation of the concentrations; and 

6. The algorithm used for solution is robust and 
allows construction of a pseudo-inverse matrix 
from basis functions that represent mixtures of 
components and that may even contain some 
noise and other background signals. 

Another embodiment 31, shown in Figure 4, 
uses the same collection optics 17 and grating 
spectrometer 27 as described above, but replaces 
the waveguide (13 in Figure 1) with a small sample 
cell 33 positioned at the focal point of the first 
parabolic mirror 17. A laser beam 34 passes 
through the sample cell 33 in a direction normal to 
the plane of the spectrometer 25 and collection 
optics 17, with the polarization of the laser beam 
34 being oriented normal to a line that connects 



30 scribe each gaseous component. 

pseudo-inverse matrices can be determined for 
these new multiple slit basis functions and used to 
calculate concentrations of the sample gases 
present. One advantage of this technique is that a 
35 plurality of slits admits more light to the spectrom- 
eter, thereby increasing the optical efficiency of the 
system. The maximum number of slits that can be 
used is determined by the resolution of the spec- 
trometer and by the patterns of the spectral lines. 
40 Slits of different widths can be combined, and the 
spacing between the slits can be variable. 

Another embodiment, illustrated in Figure 6, 
replaces the grating spectrometer 25 with a vari- 
able wavelength filter ("VWF") 51 placed over or 
45 adjacent to a photosensor array 53. The collection 
optics 17 would direct the Raman scattered light 
onto the VWF and photosensor array combination 
for wavelength dispersal and sensing. The VWF is 
disclosed in another patent application, U.S. S.N. 
50 07/ 480,172, filed on an invention by Kent Vincent 
and assigned to the assignee of this patent applica- 
tion. The VWF 51 is constructed so that the central 
wavelength of the narrow pass band of the filter 
varies linearly or strictly monotonically along the 
55 length of the VWF. Each of the elements of the 
photosensor array detects only a narrow band of 
wavelengths, determined by the position that the 
transmitted or reflected light issues from the VWF 
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51. Such a VWF can be made today with band 
pass intervals of approximately 0.5-1.0 percent of 
the central wavelength. A VWF 51 may include a 
narrow band wavelength filter 51a and a broad 
band wavelength filter 51b, positioned in tandem as 
shown in Figure 4. The group of light wavelengths 
issued (transmitted or reflected) by the narrow 
band filter 51a at any point on this filter is quite 
narrow, as narrow as 2 nm, and the central fre- 
quency for this pass band increases monotonically, 
either continuously or discontinuously, as the posi- 
tion of the point is displaced in a selected coordi- 
nate direction (x) along the filter. The wide band 
filter 51b has a somewhat wider wavelength pass 
band at any point on this filter. The central fre- 
quency of this wide band filter approximately co- 
incides with the central frequency of the narrow 
band filter 51a at any corresponding coordinate x 
for these two filters. The pass band of the wide 
band filter cuts off all light wavelengths for which 
the narrow band filter 51a would manifest unwanted 
side band emission so that the combination of the 
narrow band light filter 51a and the wide band light 
filter 51b produces a narrow wavelength band of 
light with no side band emissions; and the central 
frequency of this narrow band increases (or de- 
creases) monotonically with increase of the coordi- 
nate x. 

In the foregoing embodiments, the sample cell 
or waveguide tube 13 shown in Figure 1 has been 
illustrated as a long cylindrical waveguide tube, 
coated on its interior or exterior surface with a 
highly reflecting material that reflects the Raman- 
scattered light and directs the scattered light along 
the waveguide. In another preferred embodiment 
61, shown in Figure 7, a cylindrical tube 63 is again 
coated on its interior or exterior surface with a 
highly reflecting material 65, but a small aperture 
67 is formed in this material at a position near the 
center of the length of the tube 63. Raman-scat- 
tered light is emitted from the interior of the 
waveguide tube 63, and this radiation is received 
and analyzed by collection optics, a spectrometer 
and a light sensor array as illustrated in Figure 1. 
Because the remainder of the light-receiving sur- 
faces of the waveguide tube 63, excepting the 
aperture 67, are highly reflecting at the scattered 
wavelength \ = x R , the waveguide arrangement 
shown in Figure 7 will repeatedly scatter light emit- 
ted on or adjacent to the longitudinal axis LL' back 
into the region near the axis LL' after reflection 
from the tube walls. This always occurs when fight 
scatters from an interaction region that is coinci- 
dent with a narrow light beam travelling along the 
axis in a reflecting tube. Any region near the tube 
axis contains multiple images of light scattering 
regions positioned along or adjacent to the tube 
axis. A light beam collected from this region 
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through the aperture 67 has a greater intensity than 
tight emitted from a point along the tube axis, in 
the absence of a reflecting waveguide or coating 
65. Light that is ultimately emitted through the 
5 aperture 67 contains: (1) the direct image of light 
scattered from the region adjacent to and along the 
tube axis; (2) reflections of these direct images 
from the tube wall opposite the aperture 67; and (3) 
light from images of multiple reflections along the 
w entire length of the tube axis. The imaging prop- 
erties of the waveguide 63 are improved if the 
maximum diameter of the light beam in travelling in 
the waveguide is made as small as possible, rela- 
tive to the inside diameter of the waveguide tube. 
75 The aperture 67 in Figure 7 is formed by 

removing, or by not initially applying, the reflective 
coating 65 in the aperture region. The aperture 
itself can have any shape, although circular or 
elliptical shapes may be preferable for subsequent 
20 collection by circular collection optics of the light 
emitted from the aperture. The shape and size of 
the aperture will preferably be determined by the 
requirements of the external optical system. If col- 
lection optics are used, the aperture diameter 
^5 d aperture might be chosen to satisfy the relation 

where d tU be is the tube 
inside diameter. For improved performance, the 
aperture 67 is placed at or adjacent to the center of 
the tube axis, and the length L of the tube should 
30 satisfy L7d tube > 20. 

Figure 8 illustrates use of the tube 63, with the 
small aperture 67 therein, in this alternative pre- 
ferred embodiment of the invention. The light beam 
69 that issues from the aperture 69 is collected by 
35 collection optics 70 and focused on the entrance 
slit 23 of the spectrometer 25 that is constructed 
and operated as shown in Figure 1. Processing of 
the light beam 69 after the beam has passed 
through the spectrometer entrance slit 23 occurs as 
40 discussed in connection with the spectrometer 25 
shown in Figure 1 . 

In another embodiment 71, shown in Figure 9, 
the small single aperture 67 in Figure 7 is replaced 
by a radial band aperture 77 on a waveguide tube 
45 73 from which a Raman-scattered light beam 79 is 
emitted. Collection optics preferably surrounds the 
tube 73 and collects substantially all the light beam 
79 emitted from the radial band aperture 77. Pro- 
cessing of the light beam 79 is similar to process- 
so ing of the light beam 69 in Figure 8. 

In another embodiment 81, shown in Figure 10, 
the small aperture 67 in Figure 7 is replaced by a 
narrow, longitudinally oriented aperture 87 on a 
waveguide tube 83 from which a Raman-scattered 
55 light beam 89 is emitted. Collection optics prefer- 
ably surrounds the tube 83 and collects substan- 
tially all the light beam 89 emitted from the longitu- 
dinally oriented aperture 87. Subsequent process- 
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ing of the light beam 89 is similar to processing of 
the light beam 69 in Figure 8. 

Claims 

5 

1. Apparatus for determining the concentration of 
each of a set of N predetermined gases (N ^ 
1) that may be present in the airway of a 
patient in a single measurement, the apparatus 
being characterized by: w 

an incident light beam (15) of a selected 
wavelength X 0) chosen to promote Raman 
scattering in at least one of a predetermined 
group of gases that may be present in the 
patient's airway, the light beam being directed 75 
along a predetermined beam axis; 

a sample cell (13) to receive and tem- 
porarily hold gas (14) from a patient's airway, 
the cell being relatively transparent to light of 
wavelength x 0 in the direction of the light 20 
beam axis and being positioned to receive the 
light beam (15) so that gas in the sample cell 
is illuminated by the light beam; 

spectrometer means (23/25/27), having a 
beam entrance and wavelength dispersion 25 
means (27), for receiving a light beam (16) at 
the beam entrance, for directing the received 

dispersion means 



each of said predetermined gases having a 
known spectrum of Raman scattering of light, 
and for determining the unknown concen- 
trations Cj. 

2. The apparatus of claim 1, further characterized 
by said predetermined group of gases includ- 
ing N 2 , 0 2 , N 2 0, C0 2 , SF&, halothane (CF 3 - 
CHCIBr), enflurane (CHCIF-CF 2 -0-CHF 2 ), 
isoflurane (CF 3 -CHCI-0-CHF 2 ), desfiurane 
(CF 3 -CHF-0-CHF 2 ), sevoflurane (CH(CF 3 ) 2 0- 
CH 2 F), methoxylf lurane (CH 3 -0-CF 2 -CHCI 2 ), 
ethanol (C 2 H 5 OH), isopropanol (C 3 H 7 OH), 
ether, cyclopropane, acetone (0 = C(CH 3 ) 2 ) 
and water vapor. 

3. The apparatus of claim 1 , further characterized 
by said selected wavelength X = X 0 lying in the 
wavelength range 500 nm < X 0 < 850 nm. 

4. The apparatus of claim 1, wherein said spec- 
trometer means (23/25/27) is further character- 
ized by including a diffraction grating for dis- 
persal of said light beam (16) into light rays 
(28) having different wavelength components. 

5. The apparatus of claim 1, wherein said spec- 
trometer means (23/25/27) is further character- 




plurality of light ray bundles (28), each bundle 
containing light rays with a different, narrow 
interval of light wavelengths; 

light collection means (21) for collecting 
portions of the light beam (16) that undergoes 
Raman scattering by one or more of the gases 
(14) in the sample cell (13) and for directing 
these portions of the light beam to the spec- 
trometer means beam entrance; 

an array (29) of photodetectors, each of 
which receives one of the plurality of 
wavelength-dispersed light ray bundles (28), 
including a narrow range of wavelengths, from 
the spectrometer means (23/25/27) and issues 
a photodetector output signal representing the 
amount of Raman scattered light received by 
that photodetector; 

wavelength allocation means for receiving 
the photodetector output signal from each of 
two or more photodetectors (29) and for deter- 
mining and issuing as output signals the rela- 
tive or absolute amount of Raman scattered 
light received by each photodetector; and 

computation means for receiving the 
wavelength allocation means output signals, for 
expressing the values of the unknown con- 
centrations Cj (i = 1, 2, N) of said N 
predetermined gases as a system of linear 
equations in the unknown parameters Cj, with 



for dispersal of said figfi 
rays (28) having different wavelength compo- 
nents. 

6. The apparatus of claim 1 , wherein said sample 
35 cell (13) is further characterized by including a 

hollow cylindrical tube (63) having a reflective 
coating (65) on its inner circumference or on 
its outer circumference that is highly reflecting 
for light of wavelength X = X R , the tube having 

40 a tube length L and a tube inner diameter r 

with a ratio L/r of at least 20, the tube having a 
circumferential wall that is transparent to light 
of wavelength X = X R , and the reflective coating 
having a small aperture (67) removed so that 

45 incident light of wavelength X = X R is transmit- 

ted through a portion of the circumferential wall 
of said sample cell adjacent to the aperture. 

7. The apparatus of claim 1 , wherein said sample 
so cell (13) is further characterized by including a 

hollow cylindrical tube (73) having a reflective 
coating (75) on its inner circumference or on 
its outer circumference that is highly reflecting 
for light of wavelength X = X R| the tube having 
55 a tube length L and a tube inner diameter r 

with a ratio L/r of at least 20, the tube having a 
circumferential wall that is transparent to light 
of wavelength X = X R , and the reflective coating 
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having an aperture (77) in the form of a cir- 
cumferential band on the tube removed so that 
incident light of wavelength X = X R is transmit- 
ted through a portion of the circumferential wall 
of said sample cell adjacent to the aperture. s 

8. The apparatus of claim 1 , wherein said sample 
cell is further characterized by including a hol- 
low cylindrical tube (83) having a reflective 
coating (85) on its inner circumference or on w 
its outer circumference that is highly reflecting 

for light of wavelength X = X R , the tube having 
a tube length L and a tube inner diameter r 
with a ratio L/r of at least 20, the tube having a 
circumferential wall that is transparent to light 15 
of wavelength X = X R , and the reflective coat- 
ing having an aperture (87) in the form of a 
narrow longitudinally oriented band on the tube 
removed so that incident light of wavelength 
X= X R is transmitted through a portion of the 20 
circumferential wall of said sample cell adja- 
cent to the aperture. 

9. The apparatus of claim 1, wherein said com- 
putation means is further characterized by de- 25 
termining the unknown concentrations Cj by a 
singular value decomposition method in which 
values of said wavelength allocation means 
output signals serve as forcing functions for a 
linear system of equations. 30 

10. A method for determining the concentration of 
each of a set of N predetermined gases (N ^ 
1) that may be present in the airway of a 
patient in a single measurement, the method 35 
being characterized by the steps of: 

positioning a sample of the patient's air- 
way gases (14) in a sample cell (13) that has 
at least one sample cell wall that is transparent 
to incident light of wavelength near X = X 0 , 40 
where Xo is a predetermined wavelength lying 
in the range 500 nm < 0 < 850 nm; 

directing an approximately monochromatic 
light beam (15) having a wavelength X 0 
through the sample cell, (13) and allowing the 45 
light beam to undergo Raman scattering by the 
sample gas (14) contained in the sample celt; 

collecting the Raman scattered light (16) 
outside the sample cell (13) and dispersing 
this light into a plurality of wavelength-dis- 50 
persed light ray bundles (28); 

providing an array (29) of photodetectors, 
each of which receives one of the plurality of 
wavelength-dispersed light ray bundles (28) 
and issues a photodetector output signal repre- 55 
senting the amount of Raman scattered light 
received by that photodetector; 

determining and issuing as photodetector 



11. 



array output signals the relative or absolute 
amount of Raman scattered light received by 
each photodetector; and 

expressing the values of the unknown con- 
centrations Cj (i = 1, 2, N) of said N 
predetermined gases as a system of linear 
equations in the unknown parameters Cj, with 
each of said predetermined gases having a 
known spectrum of Raman scattering of light, 
and determining the unknown concentrations Cj 
by a singular value decomposition method in 
which values of said photodetector array out- 
put signals serve as forcing functions for the 
linear system of equations. 

The method of claim 10, further characterized 
by including the steps of: providing said sam- 
ple cell (13) with a circumferential wall (63, 73, 
83) having a reflective coating (65, 75, 85) 
thereon that is highly reflective for light of 
wavelength X = X R ; and 

providing a small aperture (67, 77, 87) in 
the reflective coating (65, 75, 85) so that light 
of wavelength X = X R that is incident on the 
aperture is transmitted through the aperture 
and through the circumferential wall. 
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